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ABSTRACT 

We first review the properties of the conventional r-functions of the KP and 
Toda-lattice hierarchies. A straightforward generalization is then discussed. It 
corresponds to passing from differential to finite-difference equations; it does not 
involve however the concept of operator-valued r-function nor the one associated 
with non-Cartanian (level k ^ 1) algebras. The present study could be useful to 
understand better g-free fields and their relation to ordinary free fields. 
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1 Introduction 



The T-functions originally associated with integrable Toda-lattice hierarchies, have been iden- 
tified as sections of determinant bundles over the universal moduli space and also seen 
to admit a group-theoretical interpretation [[|. They are now attracting attention again be- 
cause of their appearance as non-perturbative partition functions in quantum field/string 
theory ||. This is giving a new stimulus to the search of a broader interpretation of these 
r-function which would go beyond free-fermion models and commuting differential flows ||. 
The two most obvious directions of generalization consist in moving from free fermions to 
generic conformal models on one hand and from continuous 2-dimensional spectral surfaces 
to their quantum and/or discrete counterparts on the other hand. In the group-theoretical 
approach this corresponds to going beyond the level k = 1 Kac-Moody algebras and from 
"classical" to quantum groups respectively. 

The purpose of this note is to discuss briefly some preliminary material, concerning the 
"quantum T-functions" that follow from the second direction of generalization. We concen- 
trate here on a very particular aspect of the corresponding theory, i.e. on a possible relation 
of these "quantum r-function" to the conventional "classical" Toda-lattice r-function inter- 
preting the former one as a c-numbers. We proceed by rephrasing results obtained in refs. |5[], 
where solutions to finite-difference analogues of Hirota equations are given. These solutions 
are found for the "first" time-variable. We present here one particular way to introduce the 
other "times". In this construction the c-number "quantum r-function" is obtained from 
the "classical" one, merely through a change of variable which appears to be a kind (albeit 
special) of Miwa transformation. Within this framework there is no need to go beyond the 
standard Segal- Wilson infinite Grassmannian and the standard free fermions (no g-fermions 
are required). Relations to generic operator- valued objects as well as to other aspects of the 
theory of quantum r-functions will be discussed elsewhere. 

2 Conventional r-function 

We begin by reminding the reader of the definition of the conventional r-function for the Toda- 
lattice and KP hierarchies. There are actualy several more or less equivalent definitions. Any 
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generalization should prove natural irrespective of which definition is adopted. For the sake 
of completeness and future use we give below the whole list. In this note however, we will 
deal mostly with only one of these definitions. 

The conventional r-function can be characterized by any of the following properties: 

1) It satisfies a bilinear Hirota equation, which for the KP hierarchy looks like 



j> ^V + {z\t} ® V-{z\t'}r{t} ® r{t'} = 0, 
V ± {z\t} = exp f±E 

\ n>0 



1 d_\ 

nz n dt n I ' 



2) It can be written correlator of free fermions 



r N {t,t} = (N\e H Wge B &\N), (2) 



where 



H{t} = E tnJ+n, H{t\ = ^2 t n J-n, 
n>0 n>0 

J ( z )= E JnZ^ 1 = 1p{z)lp{z) = 8(f){z), Q = exp I ^ Amn^n^n 



(3) 



iez iez 
ip n \N) = 0, for n > N, ^ n \N) = for n < N . 

3) In Miwa coordinates, the KP r-function can be represented in the determinant form: 

t{(} _ det*^ (4) 

where 

tn = 4^ + -E K\ tn = *f + - E AS, (5) 

n a n ~ 

and A(A) = Y[ a> p{\ a — -\g). The analogous representation of the generic Toda-lattice r- 
function is a little more sophisticated, see refs. |S], |S| for details. Sometimes it is convenient 
to interpret the A a 's as eigenvalues of some matrix A |5|. ||. 

The Miwa coordinates can be used to introduce a sort of Fourier transform with respect 
to the variables {t,t}, 

T N {t,t} = T{p(\),p(\)}, (6) 
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with 

tn = t^ + -[ p(X)X- n , i n = i^ + - [_p(\)X n . (7) 



n Jdx n Jd\ 

If p(X) = J2aPa${^ — A a ) and all non- vanishing p a = 1, we recover the parameterization 



In terms of f(p a ,p a ) the Hirota equations become finite-difference equations |L0] : 

(A a - h)T N (p a + l,Pb + l;Pc)tN{Pa,Pb,Pc + 1) + 

(8) 

+ ( cyclic permutations of a, b, c) = for any three A a , A c ; 
T N (p+ 1,p)t n (p,p+ 1) - f N (p + l,p+ 1)t n (p,p) = 

'-r N +i(p+ l,P+ 1)tat-i(p,p) for any two A , A fe ; (9) 



A, 
A„ 



These equations are an immediate consequence of (]!]). Actually, there is an infinite set of 
equations, but it is sufficient to consider only (||) in order to fix the KP flows associated 
with positive times, the conjugate equation is necessary to fix the KP flows w.r.t. negative 
times. Moreover, it is enough to consider only (|]) in order to obtain the complete Toda-lattice 
hierarchy. 

4) Most relevant to our discussion is another determinant formula for the Toda-lattice 
r-function §,§: 

r N {t,t}= det H id {t,t}, (10) 

— oo<i,j<N 

where the matrix Hij is characterized by the property, 

k 



±H -H - ± 
dt k ij i+k ' j \dt 1 



These equations can be solved by a sort of Fourier transform, and their general solution is 
provided by an arbitrary measure function \i(z,z): 

H = f (f[ J"'") »(z,z) (f[ e^ k ) , T lm = f z l pi{z,z)z m , (12) 

Jdzdz \ k=1 J \ k=Q J Jdzdz 

which implies that 

Hij{t, t} = H Pi-i{t}T lm P m ^{t}. (13) 



Ti m is already independent of t and t (in fact, T/ m describes the adjoint action of the element 
of the Grassmannian Q on fermionic modes). The Schur polynomials P m {t} are defined by 



oo 



JJe^ "= J2P m {t}z m , (14) 

fc=0 m=0 

so that 

^-P m {t} = P m - k {t} = f p m { t }. (15) 

The object, defined by the r.h.s. of eq. fllOD , satisfies the set of Hirota equations and thus 
is a Toda- lattice r-function as a corollary of flTT| ) . However, the origin of the Hirota equations 
from this point of view is somewhat more general than eq. (|il"T). 



3 The origin of the bilinear equation 

The Hirota equations or bilinear identities have different interpretations according to the 
different constructions which are used to interpret the r-function: 

1) From the point of view of integrable hierarchies these just appear as a compact way to 
write the entire hierarchy in terms of a single identity ([!]). 

2) In the formalism of free fermions the simplest origin of the bilinear identites is the 
formula for the quadratic Casimir operator of the corresponding version of GL(oo) ® GL(oo), 
X^L_oo ^n^^-ni which, together with the embedding of GL(oo) into the universal enveloping 
of U(l) k _ v straightforwardly implies ([!]). 

3) Eqs. (|8|), (|S|) for vanishing multiplicities p a ,pb,p c ,P,P are immediate consequences of 
the determinant formula (|4|). Eqs. (§), @ for arbitrary integer p's can also be obtained in 
this way, if p a > 1 is interpreted as the number (multiplicity) of coincident "eigenvalues" A a . 
If Xb is smoothly moved along the spectral surface and reaches the position of A a , then p a 
changes abruptly: p a — > p a + pb- Within this context all vanishing p a can be considered as 
associated with A a = oo, and the change of p a by one p a — > p a + 1 is interpreted as "bringing" 
one extra "eigenvalue" A a from infinity. 

4) The most general bilinear identities can be deduced by the following procedure. Intro- 
duce two operators, M and M, acting on the infinite-dimensional matrix Hif 

M : Hh — //, .;.,: M : //, ; — //,. ; , ; . (16) 
4 



Select some n x n submatrix of H^ and consider its determinant, 

T( n *)[a,a\P,j3] = det H^, n = ft — a = ft — a. ^7) 



n<i</3 
u<j<(3 



Then M and M act on 77 n ) by shifting the a's and /3's: 

M: a — ► a + 1, /3 -> /3 + 1, a -> a, ft ^ ft; 
M : a ^ a, /? -»• /3, a -»• a + 1, j3-*j3 + l, 
and the following bilinear relation is true:[] 

(/ MA* - .1/ M) r {n) [a,a\ft,ft}r in) [a,a\ft,ft} = 
= T( n+1 )[a, a\ft + l,ft+ l]r( ft _i)[a + l,a + l\ft,ft). 



(21) 



This formula is an identity, provided the action of M and M is defined according to 



(16). With the same definition any two elements of the matrix can be related through 



H i:j = M ; ' \I J 3 E Vjly or 



H tj = M'm'h (22) 
(where H = H 00 ). Thus 

r (n) [a, a\ft, ft] = det (m'm'h) . ( 23 ) 



a<j<0 



One can now reformulate the statement and say that (23) is always a solution to fl2"Tp 



without any references to the matrix H^ and the rules (0). This is the general statement 

about the determinant solution to the bilinear equations which we will need in this paper. 
1 It is a particular (p — 2) case of the general identity for the minors of any matrix, 



Hi. : II, ...,. r .... — , / j ( ) H h 



, ^ . ip-i\jp(.i)-jp(p-i)^rj PW ) (19) 

y ' P 



where the sum on the r.h.s. is over all permutations of the p indices and if^...* ^...j denotes the determinant 
(minor) of the matrix, which is obtained from Hij by removing the rows i\ . . . i p and the columns j\ . . . j p . 



Using the fact that (H 1 )y = H^j/H , this identity can be rewritten as 

HH^ ip \j 1 j p — i— J ) (— ) ^ip' (1) ...ip, (p ^ 1) \jp(i)...jp( p - 1 )^ip, [p) \jp( p )- (20) 

' p,P' 

Taking now p — 2 and i\ — a, i^ — P, ji — ot, 32 — ft one obtains (pi]). 
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Let us note, however, that this solution has the form of the determinant of a finite matrix. 
Such objects are usually interpreted as r-functions of forced hierarchies Generic r- 

functions can be obtained from these through a limiting procedure a, a —>■ oo: 

t n = lim t m \N -n,N- n\N, N\. (24) 

We will use in the following both determinant representations since the infinite matrix one 
is natural from the viewpoint of the fermionic realization, while solution (|23| ) is convenient to 
discuss algebraic properties. 

It might be good at this point to explain, how the conventional differential Hirota equation 
can be deduced from this general construction. It arises from a particular choice of operators 
M and M. Namely, consider a function F(X, X) and define: 

MF(X, X) = F(qX, X), MF(X, X) = F(X, qX) (multiplicative action) (25) 

or 

Mf(x, x) = f{x + e, x), Mf(x, x) = f(x, x + e) (additive action). (26) 

One can also introduce an algebraic version of the group action, specified by the operator M. 
It is defined in terms of the "derivative" operator^] 

D = a(M-I), (27) 

and similarly for D. In the additive case a simply acts as multiplication by a c-number, 
a = e -1 . In the multiplicative case a is usually defined to act as multiplication by (^zW > 
(analogously, a is defined to act as multiplication by the factor (=zWy ) and has a non-trivial 
commutation relation with M: 

aM = qMa. (28) 
As a corollary of the above definitions, we have: 

n-l 

D n = [a(M -/)]" = a" J] (Mq- k - f) = a n M n q~ n{ - n ~ V>/2 (M- 1 ; q) . (29) 



fc=0 



2 Often a more symmetric definition D = uiM 1 ! 2 —M 1 / 2 ) is used. It gives expressions in a more symmetric 
form like the analogue of (||), A(£>) = D ® M 1 ' 2 + M -1 / 2 ® D 

lengthy. Since this is not important for our presentation, we stick here to the simpler definitions 



, but makes formulas somewhat more 

diag 



6 



When acting on a product of two functions F(X)G(X) (i.e. on the "diagonal" of the tensor 
product) M is characterized by the following comultiplication property: 



diag 



i.e. 



A(M) = M®M 
M[F{X)G{X)\ = F{qX)G{qX) = [MF{X)][MG{X)}. 
Since a acts on the diagonal as multiplication, 



(30) 



A(*)ldtag = *. 



(31) 



we conclude that 



A(A)L,= A(a<jfr-i))| =»A(A-/) 



diag 

a (M ®M — I®I) 



diag 



= (Dig) I + I®D + a- 1 D®b 

diag V 

(b®M + i®b 



diag 



(32) 



diag 



Since a 1 = (g— 1)X vanishes at q — 1, we obtain in this "classical" case the ordinary Leibnitz 
rule A(D) = D®I + I®D 

diag 

The bilinear identity ([?!]) can be, of course, rewritten in terms of the operators D and D 
instead of M and M. Indeed, using the action of the operators M and M (^5|), (|26|) , one can 
write: 

T {n) (qX,qX)r (n) (X,X)-r in) (qX,X)r (n) (X,qX) = r (n _ 1} (X, X)r {n+1) (qX, qX) (33) 
in the multiplicative variables, or 



T (n) (x + e,x + e)r {n) (x,x) - r {n) (x + e, x)r {n) (x, x + e) 
= r (n _i) (x, x)r (n+1) (x + e, x + e) 



(34) 



in the additive variables. It can be equally rewritten in a form similar to (|21"D: 



or, using 



/ ® MM .U M) r (n) r (n) = / ® MMr( n _i)T( n+ i), 



(l®bb-D®D) r (n) r (n) = / ® CTMaMr (ri _i ) r {n+1) . 



(35) 



(36) 
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More important, the same is true for ri n ) as given in ( p3[ ) (if it is considered as the r- function 
of a forced hierarchy - see above). Indeed, 

t$ = det WMH = det(I + a^D) 1 ^ + a^DYH = 

0<i,j<N 

= ( q q)^-m-2)/2 & -N(N-l)/2 r N(N-l)/2 

One is actually free to choose an arbitrary normalization of the r-function. In particular, a 
natural choice, which cancels the a-factors in the r.h.s. of (|36|) (compare with (||), (|9])), is 

f# } = det D'b'H. (38) 
These normalized r-functions satisfy the equation 

(l®bb-D®D) f N T N = (/ ® MM) fjv-ifiv+i (39) 

which should be compared with (^). In the q — 1 case, formula ( |3lf ) is exactly the classical 
expression of the theory of the Toda-lattice r-functions, included as p. 4) in the list of sect. 2. 
For q 1 this result was first derived in refs. and interpreted there as the determinant 
formula for solutions of a finite-difference Hirota equation. 



4 Introduction of higher time- variables 

There is still one ingredient missing in our discussion for generic q (we consider \q\ < 1 
whenever it might lead to ambiguity). The thing is that the operators M and D are defined 
as acting on functions of a single variable X. Comparison with the q = 1 case implies that X 
plays the role of the first time- variable X = T\ (and X = 7\ while N can be considered as one 
of the zero-times). According to this, (^) is only the subset of Hirota equations, associated 
with the first time-variables. A way to introduce higher times is suggested by the relation 
(|TT| ) for the q — 1 case. A natural generalization is 

D k H = {D^H, D k H = {D^H, (40) 
where b k is a finite-difference operator w.r.t. the fc-th time 7*, 

W t ) = F( ^ rt) -y rt) . (41) 
(g fc - l)T k 



We reserve the possibility to take q k ^ q = q\. For example, the choice 



q k = q k (42) 

seems to be a reasonable alternative. 

It was a corollary of ([11]) that Hij = d l d^H = didjH acquired the form (|T3|), which is 
important for establishing the relation to free fermions. Let us now derive the analogue of 
(|l~3| ) for q ^ 1. The generic solution to eqs.(|40|) is given by an integral formula involving 
g-exponentials[} 

H = [ dzdz (f[ e qk (T k z k )) n(z,z) ( g e qk (f k z k )) . (43) 

J \k=l / \fc=0 / 



The measure /j>(z, z) is not specified by 

Introduce now the g-Schur polynomials Pjf^{T} through 



T[e qk {T k z k ) = Y.P^{T}z n . (44) 

k=0 n=0 



It follows that D k P^{T} = P^l k {T} and that 

T lm = f z l fi(z,z)z m . 



(45) 



dzdz 



The next question is that of the free-fermion representation of = detH^f. It could 



ij 



at first seem natural to look for a representation in terms of g-fermions, g-oscillators etc. 
However it turns out that the most straightforward option is to stay with the ordinary free 
fermions, the same as for q — 1, and thus, with the ordinary Segal- Wilson Grassmannian. 



3 Let us remind the reader that the g-exponential e q (x) = l/E q (—x) is characterized by the following set 
of properties: 

1. D + e q (x) = e q (x); e q (x) = e q ((1 - q)x) and lim^i e q (x) = e x ; 
2 e (x) - T * E(x)-Y « Hh - 1)/2 * k 
where (a, q) k = ^(1 - a?) = 

3. e q (x) = j^-, thus E q (x) = (-i.gjoo and 6 00 (x) = YX^q^ ' 2 x k = {q,q) oE q {q 1 l 2 x)E q (q 1 i 2 X - 1 ); 

4. E q (x)E q (y) — E q (x + y) and e q (y)e q (x) = e q (x + y), provided xy = qyx; 

The first three properties explain the relevance of the g-special functions as solutions to finite difference 
equations (i.e. to various periodicity constraints), the last property indicates why these same functions occur 
in the study of non-commutative algebras and problems of quantum mechanics and quantum field theory. 
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(This phenomenon has already been observed in ref. fi~2]| .) The origin of this possibility lies 
in the fact that g-Schur polynomials and other relevant objects can be obtained from their 
q — 1 counterparts by a simple and universal change of time-variables. Indeed, 



CO 



U^(T k z k ) = Y[^ kZ , (46) 

k=l k=l 



provided the t's are expressed in terms of the T's according to 



oo oo r Pn(-\ „ \n 

E^ = E " i ""I (47) 



Thus 



Pi q) {T} = P k {t} (48) 

and 

H$f{T,f} = H l3 {t,i}, T^{T,f} = r N {t,i}, (49) 



where satisfies eqs. (11), the entire set of equivalences given in sect. 2 remaining valid. 
Because of this, rjy {T, T} can be represented as 

T { N q) {T,T} = r N {t 1 t} ( = } (N\e H ^ge R ^\N) (50) 
with some Q = exp (j2 m , n A mn ^nipn) and 

nm-TTt t {4 - 7) V 

n>0 n,A;=0 "l 1 % J 

51 

H{t} = E*n^-n = E ZTj „n\ 
n>0 n,fc=0 'H 1 Ik) 

5 Miwa transformation and discrete Hirota equations 

The bilinear identity for q ^ 1, when written in the form ([34]), is spectacularly similar to 
"the discrete Hirota equation" (^). That arises when the ordinary (q = 1) Hirota equation is 
represented in terms of the Miwa coordinates. This is perhaps a little less surprising once we 
know the relation (|47|), which, for the variable 7\, is nothing else than the Miwa transformation 
(P). One difference is that 
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is essentially a simultaneous change of multiplicities p{q~ l X\) at all the points q~ l \i, I > 
on the spectral curve, rather than at a single point Ai = ((q — l)Ti)~ . Another difference 
is that eq. (|4j) is written in terms of the additive variables x, x, while T 1; 2\ are associated 
with the multiplicative ones T x = X, 7\ = X. We shall not go into a detailed discussion of 
these differences here but will instead address another question. 

What do the higher time-variables T k correspond to in this context? Note first of all that 
there are generalizations of the discrete Hirota equation (g), corresponding to simultaneous 
variation of the multiplicities p(X) at several arbitrary points A on the spectral curve, not 



necessarily related by a "semi-periodicity" condition like (52). It is among these equations 
that we should look for counterparts of the T k variations. In fact, the relevant set of points 
A, associated with the variable T k is 

{e 2 ^ k \ k q^ l/k \a = 0,...k-l; l>0}, \ k = ((1 - q k )T k y 1/k . (53) 

(Here we find a first reason to prefer the choice q k = q k instead of q k = q since in the former 
case we have q k k = q l .) The action of D/DT k on rff can be now described as the insertion 
of the fc-fermion non-local operator 



(\ k e 2via/k ) = : J] ij(\ k e 27Tia/k ) : (54) 

a=0 

(note that only points with I = contribute): 

M Tk 4 ] ~ (N + k * (\ k e 2ma ' k ) e H Wge B & 

e H{t}g e H{t}q {\ k e 2 ^ k 



(55) 

N + k) 



Thus, we conclude that the time-variables T k are associated with certain types of Miwa 
transformations, distinguished by a particular selection of points, where the multiplicities p(X) 
are simultaneously changed. One peculiarity is the "semi-periodicity" requirement, forcing 
to introduce all the points q~ l \, < I < oo together with any A. This is a usual thing 



in g-analysis and g-free field theory (see, e.g. |12|), which can be easily formalized in terms 



of Jackson integrals. Another peculiarity is the association of the multiples of fc-th roots 
of unity with every T^-variable. This feature is, indeed, closely related to the notion of k- 
reduction of the KP hierarchy, being, in fact, somewhat "orthogonal" to that of reduction 
(/c-reduction implies that the dependence on the variable t kn is completely eliminated, while 
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Miwa transformation, associated with does instead introduce exactly this dependence and 
nothing else). 



6 Towards (/-Matrix models 

The last promising direction for further research that will be mentioned in this note concerns 
the search for "g-matrix models" . We restrict ourselves here to several preliminary remarks 
concerning the relevant eigenvalue models and the associated conformal theories |8| , leaving 
the issue of matrices and quantum groups for a more detailed presentation. 

The partition functions of the eigenvalue models are usualy represented in multiple integral 
form 



II / z i 

■ Jdzidzi 



A(z)A(z)Y[e tkZ "Y[e Iki 

i,k i,k 



(56) 



The measure fx(z, z) in this expression depends on the concrete theory and for the case of two- 



matrix model is equal to e 



-c(z-z) 2 



This integral can be easily rewritten 0, [8], [QJ as det H, 



with Hi 



X 



dzdz 



n(z, z) TTfc etkZ Ilfc etkzk which evidently satisfies fllTD and, therefore, 



Z is a Toda-lattice r-function. The simplest model arises when c — > oo, i.e. fi(z, z) = 5(z, z). It 
is called the one-matrix model and corresponds to the reduction of the Toda-lattice hierarchy 
to the Toda-chain one. The r-function of this hierarchy depends only on the difference of 
times tk — tfc and is generally described by a particular matrix Hij = H i+ j. The relevant 



matrix model integral is 



it gives rise to 



^(lttmi) 



n 



dzi 



A 2 (z)]Je tkZ i 

i , k 



H 



(lmm) 



H, 



dz 



k 



(57) 



(5? 



There is an obvious way to deform these eigenvalue integrals. One can just change all 
exponentials for g-exponentials. This is not enough, however, to clarify the connection with 
quantum groups. The following remarks makes this relation a little clearer. The thing is that 
the appearance of Van-der-Mondians in (|56| ) is crucial for the representation of such integrals 
in the form det Hy, characteristic of r-functions associated to (forced) Toda-lattice hierarchies. 
However, Van-der-Mondians are not very natural objects in the theory of quantum groups. 
Their two natural analogues are either: 
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(z) the quantum dimension which looks essentially like Yl^jiq^ Z] — q Zj Zl ) [0 and is 
unsuitable if the deformed integral still needs to be a Toda r-function, or 

(ii) an object like 



IliQZi -Zj)( 



(59) 



i>j 



which arises in the theory of g-free fields (g-affine algebras) ||14|| . Remarkably, in the course 
of integration, it is possible to change variables to replace (|59|) by A 2 (z) and and to obtain 
the appropriate determinant formula. Indeed, 



n 



dzi 



Y[(zi -Zj){qzi -Zj)J[e tkZ " 

J i>j i,k 

nli(i-g a ) i 



'l-q) N N\ 



n 



dz. 



(60) 



t k z k 



i>j 



■i,k 



with Hij defined in (|58|). 

The partition functions of the eigenvalue models can be usually represented as correlators 
in two-dimensional conformal field theory || [|. This is especially useful for the study of 
Ward identities. In particular, for the one-matrix model we have: 



(61) 



A' 



where <fr(z) is the free scalar field and J(z) = d<fi(z). 

To deform this model in the way described above, one should perform the substitution of 
times ([5l"l) 

^ _ ^ __ (1-aY „ 

(62) 



Note now that the square of the Van-der-Mondian in fl57|) arises from the operator product 
expansion of the s/(2) fc=1 -currents in (|6T|): 



n^ w =n^(*)~n(*-*i) 2 - 

i i i>j 

If these currents are replaced by those of the sl q (2) k=1 algebra fl4]], we get 



(63) 



II J q( Z i) ~ Ii( Z i ~ Z l)(<l Z i ~ Z i)- 



(64) 



We thus obtain an argument in favor of the integrable system (0). The l.h.s. of (|60D is 
reproduced in this fashion only for a particular choice of g-free fields. (It depends on the 
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precise choice of factors in front of the negative and positive harmonics of the field <f>(z); this 
can be done in different ways, without changing rules like (|5~4]).) 

A very important ingredient in the theory of eigenvalue models is the fact that one has 
explicit expressions for the Ward identities which usually take the simple and recognizable 
form of Virasoro and W-constraints.f] In our g-matrix models, the Ward identities should 
presumably obey some kind of q- Virasoro and q-W constraints. We shall not go into any 
details here, we shall merely note one important identity, that arises if the choice 



Qk = q 



; (65) 



is made. This choice respects the fact that the argument of the ^-exponential, while linear in 
Tfc, involves the fc-th power of z. Because of this 

D . fT k , e qk (T k z k q k ) - e qk (T k z k ) 



Dz qKK ' (q-1 



(66) 

T k-^re q u(T k z k ). 



q—1 DTy Cq 



■k 

This does not immediately lead to some simple formulas for the constraints, because of the 



twisted nature of the g-Leibnitz rule: in actual calculations, for models like fl5"7f) , every D/Dz a 
will be acting on the product of e q k(T k z k ) with different fc's, while every D/DT k will be acting 
on the product of e qk {T k z k ) with different a's. Integration by parts is also somewhat more 
complicated than in the q = 1 case. To make this integration possible, one should certainly 
use the adequate notion of integration, inverse to that of finite-differentiation. It is essentially 
provided by the Jackson integral, 

! F{X) = {l-q)Y J F{q l )q l . (67) 

Jq l 

Usually the sum on the r.h.s. is defined over all integers. Then, the integral of a total 
derivative is given by the boundary terms, J DF(X) = F(oo) — F(0). In order to get q- 
Virasoro constraints, it would be useful to have the relation J q DF(X) = 0. For this purpose 



it will not be enough to simply use the integral (|67|), the reason being that it tends, in 



4 The fact that the full set of Ward-identites (constraints) forms this kind of algebras is very general, at 
least, as general as the existence of some free-field representation of partition functions, in turn equivalent to 
their identification as r-functions. It is remarkable that in the simplest matrix models these algebras arise 
naturally in some simple and well known representations. 
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the classical limit, to an integral over the semi-axis (0, oo) and that there is no need for a 
reasonable function F(X) to vanish at zero. 

There is however another integral which in the classical limit goes to an integral over the 
entire real axis, namely, 



A simple example of a function, which decreases rapidly at both ends of the integration 
domain (i.e. as X -> ±00) is e q (-X 2 ) = Y\™ =0 (1 + q k X 2 )' 1 . 
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